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ABSTRACT 



A star cluster in a galactic nucleus sinks toward the galactic center due to dynamical 
friction. As it spirals inward, the cluster loses mass due to stellar evolution, relaxation driven 
evaporation, and tidal stripping, eventually dissolving in the galactic tidal field. We model the 
inspiral of dense young star clusters near the center of our Galaxy to study the extent of the 
region of parameter space in which the cluster can reach the inner parsec of the Galaxy within 
a few million years. Since we neglect changes in cluster structure due to internal evolution, the 
present study is most applicable to star clusters less than about one initial relaxation time old. 
We find that only star clusters with initial masses 10^ Mq can reach the Galactic center from 
an initial distance of ^ 60 pc within one initial relaxation time or a few million years, whichever 
is smaller. 

Subject headings: stellar dynamics — methods: analytical — Galaxy: Bulge — 
Galaxy: Nucleus — Galaxy: Stellar content — open clusters and associations: 
individual IRS 16 — 



1. Introduction 

The innermost 100 pc of the Milky Way Galaxy contains a number of intriguing objects. These 
include the central ^ 2 — 3 x 10^ black hole (Genzel et al. 2000; Ghez et al. 2000), a cluster containing 
at least 15 massive young stars (Tabmley & Rieke 1993; Krabbe et al. 1995), a much larger population of 
older stars (Alexander 1999), and at least two young dense star clusters — the Arches and the Quintuplet 
systems (Nagata et al. 1995; Nagata et al. 1990; Okuda et al. 1990). 

Krabbe et al. (1995) found ^ 15 bright He I emission line stars in the Galactic center. They are part of 
the co-moving 7-8 Myr old > IO^Mq association known as IRS 16 (Tamblyn & Rieke 1993; Krabbe 1995), 
and are accompanied by many less luminous stars of spectral types O and B (Genzel et al. 2000). Detailed 
spectroscopic analysis of the Galactic center region (Najarro et al. 1994) indicate that these emission-line 
stars are evolved, with a high surface ratio of helium to hydrogen riHc/'^H = 1 to 1.67. Allen et al. (1990) 
classify them as Ofpe/WN9 stars, while Najarro et al. (1997) identify them as 60-100 M0 Luminous Blue 
Variables (LBVs), the late evolutionary stages of very massive stars (Langer et al. 1994). Depending on 
the interpretation of the data, the age of IRS 16 therefore lies in the range 3-7 Myr, the lower figure 
corresponding to the LBV identification. 
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One possible explanation for these stars is a recent ~ 10^ M© starburst (Krabbe et al. 1995). However, 
this model is problematic, as the formation of stars within a parsec of the Galactic center is difficult; the 
Galactic tidal field is sufficient to unbind gas clouds with densities ^ 10'' cm^"^ (Giisten & Downes 1980). 
Gerhard (2001) has proposed that a million-solar-mass star cluster formed at a distance of ^ 30 pc from 
the Galactic center could have reached the Galactic center via dynamical friction before being disrupted by 
the Galactic tidal field or by internal dynamical evolution. This qualitative argument solves the problem of 
the presence of young, very massive stars in the Galactic center. Gerhard's dynamical friction time scale 
assumed that the stellar density in the vicinity of the Galactic center is described by an isothermal sphere; 
in addition, he ignored stellar mass loss and the internal dynamical evolution of the cluster. In this paper 
we present a more quantitative approach to the problem. 

This is the first in a series of papers in which we consider the time scale on which a star cluster 
sinks to the Galactic center and is disrupted by the Galactic tidal field. In the semi-analytic calculations 
presented here, we study the inspiral of three quite different cluster models. We begin with the simplifying 
approximation that the inspiraling object has constant mass. Later, we relax that assumption and allow the 
cluster to lose mass, first by tidal stripping, then also by stellar evolution and relaxation. For definiteness, 
and for purposes of illustration, we adopt a simple analytic prescription for mass loss from the evolving 
cluster, and investigate its consequences. In a future paper we will incorporate more realistic treatments of 
cluster dynamics. 

The organization of this paper is as follows. In §2 we first consider the orbital decay of a unevolving 
point mass. In §3 we expand our study to include clusters of nonzero radii, allowing their masses to evolve 
in time as material is stripped by the Galactic tidal field. The introduction of physical parameters into 
our models then allows us to incorporate simple treatments stellar mass loss and evaporation within our 
simple model. In §4 we apply the model to star clusters near the Galactic center, to determine the region 
of parameter space in which clusters can transport a considerable fraction of their initial mass to within a 
few parsecs of the Galactic center before disruption. We discuss our results and conclude in §5. 



We begin our study with the simplifying assumption that the mass of the inspiraling object is constant. 
This idealization may be appropriate for a single massive black hole or a very compact star cluster which 
is much smaller than its Jacobi radius, the limiting radius of a cluster in the tidal field of the Galaxy. In 
the latter case, however, for the constant-mass approximation to hold, internal dynamical evolution of the 
cluster should also be negligible on the time scale on which the cluster sinks to the Galactic center. In 
practice, especially for small clusters, this will not be the case, as we discuss in §5. 



We characterize the mass M within a sphere with radius R centered on the Galactic center as a power 

law: 



2. Inspiral with constant mass 



2.1. Dynamical friction 



M{R) = AR" , 

where A and a are constants, with 1 < a < 2 of interest here. The density at distance R then is 



(1) 



>a—3 



(2) 
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and we can write down expressions for the orbital acceleration at distance R from the Galactic center 

a{R) = GAR°'-'^, (3) 

the potential 

CA 

m = — 7i?"-\ (4) 

a — L 

the circular velocity 

Vc^Ci?) = GAR''-^, (5) 

and the total energy of a circular orbit 

E,{R) = IGAR-^-^ (^) • (6) 

The object's acceleration due to dynamical friction is (Binney & Tremaine 1987, p. 425) 

aj = -47rlnAGVrn^X- (7) 

Here, m is the mass of the object, Vc is its velocity vector (in a circular orbit around the Galactic center), 
InA ~ lii(r)/i? ~ 5 is the Coulomb logarithm (where (r) is the object's characteristic radius — roughly the 
half-mass radius in the case of a cluster), G is the gravitational constant, and 

X = erf(X)-^e-^^ (8) 

where X = Vc/V^a and cr^(ii) is the local one-dimensional velocity dispersion, assumed isotropic. 
Substitution of Eqs. (1), (2) and (5) into Eq. (7) results in 

«/ = = "XlnA^, (9) 

from which we note that 

For a = 1.2, we obtain X ~ 0.89 (see Appendix A), and hence x ~ 0.34. With InA = 5 we find 
af/a{R) ~ 2m/ M. 



2.2. Orbital decay 

We can now derive the inspiral time scale for a star cluster with constant mass m in a power-law 
density profile given by Eq. (1). The time derivative of Eq. (6) is 

= -,l„AG=^, (11) 

Vc 

where the second equation expresses the work done by dynamical friction (Eq. 9). Hereafter, R should be 
interpreted as R{t), the distance from the cluster in question to the Galactic center. Substitution of Eqs. (2) 
and (5) leads to 

^ = -axlnAG3/2Ai/2TOi?(«-5)/2, (12) 
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whence 

^ = , (13) 

where 

, = 2mlnA-^(§)^^ (14) 
Solving Eq. (13) with R{t) = Rq at time t = results in 

2/(a+3) 



R{t) = Ro 



2jj(a+3)/2^_ 



(15) 



Setting R = a,t t = tdf and substituting Eq. (1) yields 



a+1 1 fMoV^^Rl^'^ 



~ a(a + 3) xln A \g) m ' ^^^^ 

where Mq — M{Rq). In terms of the orbital period of a circular orbit around the Galactic center at distance 
Rq, To = 27r (GMo/i?g)"^^^ Eq. (16) becomes 

idf ^ a+l 1 Mq 

To 27ra(a + 3) xln A m ' ^ 
For a = 1.2, Mo/m = 10^, and InA = 5, we find tdi ^ 40To. 



3. Clusters with variable mass 

We now consider the possibility that the mass of the cluster varies with time: m = m{t). Most mass 
loss from the cluster is the result of tidal stripping as the cluster sinks toward the Galactic center. We begin 
by determining the Jacobi (tidal) radius rj of the cluster in the tidal field of the Galaxy. 



3.1. Mass of a tidally limited cluster 

The differential acceleration at distance rj from the center of the cluster is obtained from Eq. (3): 

Aatide « (a - 2)GAi?"-3rj , (18) 
or, relative to the internal cluster acceleration at rj, 

l^«-el-(,_2)fil)f!2)\ (19) 



a J \ TO J J \ R 

Here aj = Gmj/rj^ and toj, the cluster mass within radius rj (still to be determined), will henceforth be 
identified as the cluster mass. Setting | Aatide | = cij, wc find 
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This may be conveniently (and conventionally) expressed in terms of average densities pj = 3mj /47rr and 
Pg = 3M/47ri?3, as 

Pj = (2-a)pG. (21) 

To proceed further, wc must make a connection between toj and rj. Two particularly simple cluster 
density profiles lend themselves easily to analytic development: 

1. A homogeneous sphere of mass mo, radius b, and uniform density 

2. A Plummer (1911) model of mass mo and scale radius b, with 

p{r)=po{l+ryby'^\ (23) 

where po is again given by Eq. (22). 

Note that, in each case we assume fi,xed parameters mo and b — that is, we neglect structural changes in the 
cluster due to dynamical evolution or stellar mass loss. This assumption greatly simplifies the calculation, 
but clearly is of questionable validity when the internal dynamical time scales are comparable to the inspiral 
time (see §5). In the next subsection we expand our model to allow for the effects of mass loss due to stellar 
evolution and escaping stars. A more complete treatment of the cluster's structural evolution will be the 
subject of a future paper. 

For the homogeneous sphere (case 1), the desired relation between mj and rj is simple: 

^^^[moir^/bf (o<6), 
[ mo (rj > b). 

No solution to Eq. (21) exists for pq > po/{2 — a), and the cluster is destroyed at Galactocentric radius 

1 l/(a-3) 

mo ' 
.(2-a)A63_ 

Outside i?min, rj >b and m = mo. Inside, rj = m = mj = 0. 
For the Plummer model (case 2), 

Pj = Po(l + raV^^)"'/' (26) 
and again, no solution exists for R < iimin- Outside iJmin, rj{R) satisfies 

/ \ 2/3 / D \ 2 / D X 2-2a/3 

: + .,7,^ . (2-.)-(i|) (t) (I) ■ ™ 

The mass of the cluster is then given by Eq. (21): 

mj(i?) = (2-a)Mo(^y (-^y . (28) 
We use this model as the basis for our discussion in the remainder of the paper. 



-^min — 



(25) 
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3.2. 



Mass loss from stellar evolution 



Many clusters dissolve so quickly that stellar evolution barely affects their mass. However, if the 
cluster survives for more than a few million years, mass loss from the most massive stars may become 
important (sec McMillan 2003 for a recent review) . Recent detailed A^-body simulations by Portegies Zwart 
et al. (2001) have quantified the expansion of a tidally-limited cluster as its mass decreases. The expansion 
drives more rapid disruption, while the mass loss slows the inspiral. 

We include stellar mass loss in our model as follows. First we rewrite Eq. (28) as mj(i?) = (m)nj(i?), 
where nj (R) is the number of stars within the Jacobi radius and (m) is the mean stellar mass, which is 
now a function of time due to stellar evolution. We assume that the mass functions of the cluster and of 
the escaping stars are identical. (Again, this is equivalent to the neglect of internal dynamical evolution.) 
We parametrize the cluster's expansion in response to stellar mass loss by 



which is equivalent to the assumption that the cluster loses mass adiabatically, as found by Portegies Zwart 
et al. (2001). 

The mean mass in the cluster can be computed from the initial mass function. For clarity we assume 
that all the mass in stars having masses above the cluster's turn-off mass is simply lost from the cluster. So 
long as the turn-off mass exceeds ^ 8M0 this assimiption is justified by the high- velocity kick imparted to 
compact objects by the supernovae in which they form, allowing them to escape from the cluster. For older 
clusters this assumption breaks down as lower- mass stars turn into white dwarfs, which do not receive high 
velocities at formation, although such clusters are not of direct interest in the present paper. Integrating 
the initial mass function, we find 



Here x is the exponent for the (assumed) power-law mass function (Salpeter: x = 2.35), mto is the cluster 
turn-off mass, and mmin is the lower mass limit. We determine the turn-off mass using fits to the stellar 
evolution models of Eggleton, Fitchet & Tout (1998). 



A tidally limited star cluster in the tidal field of the Galaxy will also lose mass due to internal relaxation 
as occasional interactions between cluster members result in velocities high enough for stars to escape the 
cluster potential. Portegies Zwart & McMillan (2002) have simulated star clusters near the Galactic center; 
they derive the following approximate expression for relaxation-driven cluster mass loss: 



6 = 60 (TO)o/(m), 



(29) 




(30) 



3.3. 



Mass loss due to relaxation 




(31) 



Here is the relaxation time at the cluster tidal radius 




(32) 



where nj is the number of stars contained within the Jacobi radius. 
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The clusters in the study of Portegies Zwart & McMillan (2002) did not spiral in to the Galactic 
center, so the relaxation time at the tidal radius remained constant over the lifetime of the cluster. In our 
case, where clusters sink toward the Galactic center, the relaxation time at the tidal radius changes with 
time. We therefore recast Eq. (31), as follows. Differentiating Eq. (31) with respect to time, identifying m{t) 
with mj and t^i with the instantaneous relaxation time at rj, and including the radial dependence of the 
relaxation time, assuming a tidally limited cluster, we obtain 

dm mo / i? 

'dt ~ ~0.29trt,o 

For the purposes of this paper, we draw a distinction between the processes of tidal stripping, in 
which stars outside the Jacobi radius are removed by the Galactic tidal field as the cluster sinks toward 
the Galactic center and the Jacobi radius shrinks, and evaporation, in which stars are driven across the 
instantaneous Jacobi radius by internal two-body relaxation. All models discussed in the following section 
include tidal stripping; models discussed in §4.2 and subsequently also include both evaporation-driven and 
stellar-evolution mass loss. 




4. Results 

Prom §2.2, the distance from the cluster to the Galactic center satisfies 

^ = -j(R) (34) 

(Eqs. 13 and 14), with R = Rq at t = 0. Transforming to dimensionless variables ^ = R/Ro and r = t/To, 
and substituting Eq. (14), we rewrite this equation in the form 

^ = i^;,lnA^r^-^^/^. (35) 
dr a + 1 Mo 

For a Plummer model, the cluster mass mj varies as a function of R and therefore ^ via Eq. (28). We solve 
Eq. (35) numerically, as it admits no simple analytic solution. For all models we adopt A = 4.25 x 10^ M© 
and a = 1.2 (Sanders & Lowinger 1972, Mezger et al. 1999). 

For simplicity we assume that x In A = 1 for the remainder of this section, unless indicated otherwise. 
A value of xlog A = 1.2 or 1.3 is probably more appropriate (Binney & Tremaine 1987; Spinnato, Fellhauer 
& Portegies Zwart 2003). The dynamical friction time scale is inversely proportional to x^^^ (see Eq. 16), 
so the effects of different choices can be easily estimated. 



4.1. Solutions without stellar evolution 

For systems without significant stellar mass loss or evaporation, the evolution may be conveniently 
parametrized by the dimensionless quantities 

/? = bo/Ro, 

H = mo/Mo. (36) 

The contours and gray shades in Fig. 1 present the dissolution time (r) of the cluster as a function of (3 and 
II. 
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Fig. 1. — Scaled cluster lifetime (contours and gray shades) as a function of the dimensionless parameters 
/3 and /u (see Eqs. 35 and 36), for models with tidal stripping but without additional mass loss by stellar 
evolution or evaporation. The parameter (3 is the ratio of the initial cluster length scale to the initial distance 
to the Galactic center. The parameter jj, is the ratio of the initial cluster mass to the mass of the Galaxy 
contained within the initial orbit. The numerical labels on the contours give the disruption time r in units of 
the cluster's initial orbital period around the Galactic center. The gray shades provide the same information 
as the contours; darker shades represent longer cluster lifetime. The dotted line indicates the values of /? and 
jjL corresponding to 6q = 0.9rj: the characteristic scale of the initial Plummer model is 90% of the cluster 
Jacobi radius. No initial solution exists for the area to the right of the curve 6o = '"j (also indicated). 

Fig. 1 shows that compact, massive clusters have the shortest lifetimes, that the lifetime decreases with 
increasing mass at fixed initial cluster radius (&o)j increases with increasing radius at fixed mass, and is 
largely independent of the radius for small radii. This last point simply means that clusters initially well 
inside their Jacobi radii (6o ^ 0.9rj) experience significant stripping only near the end of the inspiral 
process. There is no initial solution when Pq > po/(2 — a), i.e. when Jjq > rj. 

Since the stellar density diverges toward the Galactic center, no extended cluster can actually reach 
R = (although a black hole can). Fig. 2 shows the cluster's distance to the Galactic center as a function 
of time (again in units of the initial orbital period of the cluster around the Galactic center) for several 
selected values of /3 and /U. Not surprisingly, more massive clusters (larger values of fj,) spiral in more 
quickly, and physically larger clusters (larger /3) dissolve at larger distances from the Galactic center. 

The long lifetimes of clusters with bo ~ O.Qrj (see Fig. 1) and small values of n are due to the weak 

effect of dynamical friction in those cases. Since we ignore stellar mass loss and internal dynamical 
evolution (specifically, evaporation) in this simple model, such clusters survive for unrealistically long 
times. In practice, these systems will be strongly affected by stellar evolution and evaporation, as we now 
demonstrate. 



4.2. Evolution with stellar mass loss and relaxation 

By selecting the Galactic center as representative nucleus we can attach physical units to the selected 
values of fi and (3. The advantage of introducing physical parameters is that the numbers become more 
intuitive, but of course we lose the scale-free solution from previous section. Another advantage of fixing 
the scaling is that we can take stellar evolution and internal relaxation into account. Stellar mass loss (via 
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Fig. 2. — Scaled distance to the Galactic center R/Rq as a function of scaled time I/Tq, for selected 
combinations of /? and fi. The three families of models shown have /U = 0.012 (left set of curves), ^ = 0.0012 

(middle), and /i = 0.00012 (right-most curves). The dotted lines give the evolution for a constant point mass 
(/3 = 0); other curves present models with /? as indicated. The model corresponding to /x = 0.0012, /? = 0.025 
(heavy solid line) is the basis for Figs. 3-5. 



Eq. 30) and evaporation (via Eq. 33) are included by solving Eq. (35) . For most calculations we adopted 
a Salpeter initial mass function between 0.1 Mq and 100 Mq. The effect of relaxing this assumption is 

illustrated in Fig. 5 below. 

Fig. 3 shows distances to the Galactic center as functions of time of model clusters having initial 
masses of 64,000 (a) and 256,000 (b). For each selected initial distance (i?o = 2pc, 4pc, 8pc and 
16 pc) we choose a range of initial values for the cluster scale 6o = 0.2 pc, 0.4, and 0.8 pc. The choice of 
mo = 64, OOOM0, i?o = 8pc, and 60 = 0.2 pc corresponds to the "standard" model indicated in Fig. 2. 
These models were computed taking both stellar mass loss and evaporation into account. Models with 
&o = (point mass case, without stellar mass loss or evaporation) are also included for comparison; they are 
identical to the calculations presented in Fig. 2. 

Due to the extra mass-loss channels (stellar evolution and evaporation) and the resultant reduction in 

the inspiral (and hence tidal stripping) rate, the lifetimes of the clusters shown in Fig. 3 may be either longer 
or shorter than those of clusters in which stellar evolution is neglected (as in Fig. 2). This is illustrated in 
Fig. 4 for models having n = 0.00012 and /3 = 0.05, 0.025 (the "standard" model) and 0.0125. Dimensionless 
times are converted to megayears using an orbital period of 0.30 Myr, appropriate to a cluster at an initial 
distance of 8 pc from the Galactic center. The dotted curves show the evolution of the dimensionless models 
in which only tidal stripping is included. The solid, dashed and dash-3-dotted curves present the same 
models with stellar evolution and evaporation taken into account. At a distance of 8pc from the Galactic 
center, the values of /? = 0.05, 0.025, and 0.0125 correspond to 60 = 0.4 pc, 0.2 and 60 = 0.1 pc, respectively. 

Compact clusters (bo ^ 0.4 pc) are relatively luiafFccted by tidal stripping. As a result, the primary 
effect of stellar mass loss is simply to decrease the inspiral rate, increasing the cluster lifetime. However, 
in larger clusters {bo ^ 0.4 pc) the expansion caused by stellar mass loss greatly increases the stripping 
rate, significantly decreasing the lifetime despite the slower inspiral. For very low-mass, or very large, 
clusters relaxation and evaporation may dominate. However for the cases studied here the effect is almost 
negligible. This is illustrated by the first few million years of evolution of clusters with large bo, such as the 
dash-3-dotted curve in Fig. 4. The small deviation of the 60 = 0.8 pc curve from the /? = 0.050 dotted curve 
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Fig. 3. — Time-evolution of Galactocentric distance R for models with (a) m = 64, 000 and (b) 
256, 000 M0, assuming xlnA = 1. Initial cluster scales bo are 0.0 (dotted lines), b = 0.2 (solid line), 0.4 
(dashes) and 0.8 (dash-3-dotted lines). The corresponding values for fi are indicated near the start of each 
family of curves. The "standard" model marked in Fig. 2 is also indicated here (heavy solid line). For 
reference, the orbital periods of clusters at 2pc, 4pc, 8pc and 16 pc from the Galactic center are 0.09 Myr, 
0.16 Myr, 0.30 Myr and 0.56 Myr, respectively. 




Fig. 4. — Mass as a function of time for a 64,000 Mq star cluster at a distance of 8 pc from the Galactic 
center (/i = 0.00012) with various values for 60 (as indicated). The dotted curves present the corresponding 
evolution of the cluster without stellar evolution or evaporation. The values of /3 correspond to the choices 
of bo . All calculations were made using x In A = 1 . 

Fig. 5 illustrates how varying the cluster initial mass function alters the time evolution of its 
Galactocentric radius and mass. The dotted curve (Fig. 5a only) shows the constant point-mass case for 
l_i = 0.0012. The solid curve shows the evolution of the standard model with /? = 0.025, scaling times to 
megayears assuming an initial Galactocentric distance of 8pc. The dashed curves give the results when 
stellar mass loss and evaporation are taken into account, assuming lower mass limits for the initial mass 
function of 0.1 Mq, O.2M0, O.4M0 and 0.8 Mq. 

Increasing the low-mass cutoff in the mass function increases the effective stellar evolution mass-loss 
rate, and reduces the cluster lifetime. A similar effect can be achieved by increasing the power-law slope x 
of the mass function. The models are therefore degenerate in the a;-mmin plane. For example, a model with 
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Fig. 5. — Evolution of (a) the distance to the Galactic center and (b) mass for model clusters with initial 
mass mo = 64, 000 M©, initial radius &o = 0.2 pc, and initial Galactocentric distance -Ro = 8pc. The dotted 
curve in (a) shows the evolution of the model without stellar evolution or evaporation, assuming that the 
cluster is a point mass. The solid curves are the "standard" model, computed without stellar evolution or 
evaporation, with /3 = 0.025, appropriate to the choice of ho. The dashed curves include stellar evolution 
and evaporation and arc computed using a Salpeter initial mass function with different lower mass limits 
'Timin, ranging from 0.1 Mq (heavy dashed line) to 0.8 (as indicated). As before, we assume xlnA = 1. 

a Salpeter initial mass function {x = 2.35) and a low-mass cutoff at 0.2 Mq evolves almost identically to a 
model with a; = 2.13 and mmin = 0.1 or with x = 2.90 and mmm = 0.4. 



4.3. Comparison with Kim et al. (2000) 

Kim (2000) used GADGET, the tree code developed by Springel, Yoshida & White (2000), to compute 
the dynamical friction of dense star clusters near the Galactic center. In these calculations, the inner part 
of the Galaxy was represented by 2 million point particles distributed as a truncated softened power-law 
similar to our Eq. (2), except that the overall density was 2.5 times smaller than ours. The black hole in 
the Galactic center was represented as a single particle. The star cluster was modeled as a Plummer sphere 
with bo = 0.85 pc, using 10^ point particles having a total mass of mo = 10^ M©. Initially, the cluster was 
placed in a circular orbit at a distance of Rq = 30 pc from the Galactic center. The simulations ignored 
mass loss by stellar evolution and evaporation. 

The time-dependence of the cluster's Galactocentric distance, as determined by Kim (2000), is shown 
in Fig. 6. His cluster orbits become slightly eccentric during the evolution, but this seems to have little 
effect on the dynamical friction timescale. For clarity we do not show the actual results reported by Kim, 
but instead match his initial conditions, which are plotted in Fig. 6 as the rightmost dashed curve. This 
model is computed without stellar evolution or evaporation, as in Kim's simulations. Our model closely 
reproduces Kim's results when we adopt InA = 3.7, which is close to the value used by Kim. For reference, 
we also plot the evolution of a constant point-mass model (dotted curve), and a model in which stellar 
evolution and evaporation are taken into account (solid curve). To guide the eye we also plot the same 
series of runs with InA = 10. 



' 

5 10 15 

t [Myr] 

Fig. 6. — Time evolution of Galactocentric distance i? of a star cluster having an initial mass m = 10^ M0. 
The rightmost set of curves represents model 1 of Kim (2000), and have In A = 3.7; the leftmost curves have 
InA = 10. Dotted lines arc for a point-mass {0 = 0) chistcr with /i = 0.00845. Dashed lines arc for models 
with 6o = 0.85 pc, but excluding stellar evolution and evaporation, as assumed by Kim (2000). Solid curves 
include stellar mass loss and evaporation, with a Salpeter initial mass function between 0.1 and 100 M©. 
Although Kim (2000) continued his calculation for only about 9 Myr, his curve and our rightmost dashed 
curve are virtually indistinguishable. 

5. Discussion 

In this section we discuss some consequences of our semi-analytical calculations. In particular, we 
consider Gerhard's (2002) conjecture, discussed in §1, that IRS 16 and the associated young stars observed 
in the Galactic center may have been deposited there by the inspiral and disruption of a much more massive 
system. We take the two known Galactic center clusters — the Arches and Quintuplet — as templates. 

Table. 1 presents the observed parameters for the Arches and Quintuplet systems. The final columns 
give the clusters' half-mass relaxation time and the time required to reach the Galactic center, according 
to Eq. (32). (The half mass relaxation time is computed by substituting rhm for rj in that equation.) It is 
clear that neither cluster will reach the Galactic center within the next few megayears, and that both were 
probably born at roughly their present distance from the Galactic center. For these calculations we have 
again adopted xlogA = 1. 

Table 1: Observed parameters for the Arches and Quintuplet star clusters. Both lie within 35 pc (in 
projection) of the Galactic center. The first two columns give the cluster name and references, followed 
by the distance to the Galactic center, age, mass and half mass radius. The last three columns give the 
two-body relaxation time at the half-mass radius, the expected time to disruption and the inspiral time scale. 



Name 


ref. 


Rgc 


Age 


M 


Rhm 


^rt ^diss 








[pc] 


[Myr] 


[Mq] 


[pc] 


— Myr -— 


[Gyr] 


Arches 


a 


30 


2 4 


12 50 


0.2 


12 60 


0.3 1.9 


Qiiinluijlcl 


b 


35 


3 5 


10 K) 


0.5 


12 60 


5.5 9.6 



References: a) Figer et al. (1999a); b) Glass, Catchpole & Whitelock (1987); Nagata et al. (1990); Figer, 
Mclean & Morris (1999b). 
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(a) R [pc] (b) R [pc] 

Fig. 7.—- (a) Inspiral time and (b) final distance to the Galactic center for clusters with bo = 0.2 pc, as 
functions of initial galacto centric distance and cluster mass. The dotted lines correspond to bo = 0.9rj 
and bo = rj, as in Fig. 1. The solid curves indicate initial conditions where the dynamical friction inspiral 
time scale tdi is 1 Myr (left), lOMyr, and 100 Myr (right). The dashed curves in panel (a) correspond to 
^df = 0.2tih (left) and tdf — t^h (right), where t^h is the initial half mass relaxation time of the cluster, 
obtained by substituting rhm for rj in Eq. (32). The approximate locations of the Arches and Quintuplet 
clusters (see Tab. 1) are also shown. 

Fig. 7 presents, as a function of initial cluster mass and galactocentric radius, the time taken for a star 
cluster with bo = 0.2 pc to reach the Galactic center (a), and the distance from the Galactic center at which 
the cluster dissolves (b). Contours and greyscale represent inspiral time in (a) and dissolution distance in 

(b) . The cluster is deemed to have dissolved when it comes within 1 pc of the Galactic center, or when it 
has lost 99% of its initial mass. The dotted linos have the same meanings is in Fig. 1. 

The rightmost dashed line in Fig. 7 marks initial conditions for which the inspiral time scale equals the 
initial relaxation time. To the right of this curve, the cluster will experience significant internal dynamical 
evolution before disrupting. Our simple description of the cluster's internal structure is therefore unreliable 

to the right of this curve, but our expression for the evaporation rate is still valid. The left dashed curve 
corresponds to an inspiral time scale of 0.21^^^, where t^h is the initial relaxation time at the half-mass 
radius. This is roughly the core-collapse time for a system with a realistic initial mass function in which 
stellar evolution is relatively unimportant (see Portegies Zwart & McMillan 2002). Clusters with initial 
conditions to the left of the left dashed curve are thus expected to dissolve in the Galactic tidal field before 
experiencing core collapse. 

During and after core collapse (to the right of the left dashed curve in Fig. 7a) the structure of the 
cluster changes considerably, and our simple prescription for cluster disruption is unlikely to hold. We 
expect that the structural changes in these clusters will cause their dense cores to survive for longer, and 
that they will sink slightly closer to the Galactic center than indicated in Fig. 7(b) (see Gerhard 2002 and 
Portegies Zwart et al 2003 for further discussion). The change in disruption radius is not expected to be 
great, however, as the residual core masses are small and their inspiral correspondingly slow at late times. 

Fig. 7 clearly indicates that the two known nuclear star clusters, the Arches and Quintuplet, will 
not reach the Galactic center. Their survival times are determined by internal relaxation rather than by 
dynamical friction (see also Tab. 1). The inspiral time scale for these clusters is of the order of 1 Gyr, 
compared to their predicted lifetime of 100 Myr, based on A^-body simulations (Portegies Zwart et al., 2001). 
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Fig. 7(b) indicates that clusters with M ^ 20,000Mq barely evolve in Galactocentric radius on this time 
scale, but instead dissolve in situ, due to the combined effects of evaporation and stellar mass loss. 




100 




100 



Fig. 8. — Contours and greyscale map of (a) inspiral time and (b) final distance to the Galactic center for 
clusters with mo = 64, 000 M©, as functions of initial Galactocentric distance -Ro and cluster size bo- The 
dashed and dotted curves have the same meanings as in Fig. 7, and the two crosses mark the estimated 
locations of the Arches and Quintuplet clusters (see Tab. 1). 



Fig. 8 gives the time taken for a star cluster with mo = 64,000Mq to reach the Galactic center and 
the Galactocentric distance at which the cluster dissolves. In Fig. 7 the initial cluster mass m,o was varied 
at constant ^q. Now we vary bo keeping mo constant, thereby providing a second slice through the same 
parameter space as in Fig. 7. As in Fig. 7, the approximate locations of the Arches and Quintuplet clusters 
arc indicated. It is clear that small shifts in cither figure will not alter the basic conclusion that both 
clusters will dissolve at large distances from the Galactic center. 

From Figs. 7 and 8 it is clear that only massive ( ^ 10^ M©) star clusters can transport a significant 
fraction of their mass to the vicinity of the Galactic center within a few megayears. Also, even a million 
solar mass star cluster will require several tens of megayears to reach the Galactic center from an initial 
distance of ^ 30 pc. The most promising candidates to reach the central parsec of the Galaxy within 
lOMyr, but after significant mass segregation has occurred, are star clusters with masses ^ 10^ M©, born 
within about 20 pc of the Galactic center, with half mass radii of 0.2-0.4 pc. Less massive clusters, 
clusters farther from the Galactic center, or smaller (larger) clusters have greater difficulty reaching the 
Galactic center before disruption (core collapse). 

We therefore conclude that, if they originated in a massive star cluster, the stars in IRS 16 were born 
in a ^ 10^ M0 cluster at a Galactocentric distance of ^ 20 pc. The cluster deposited about 10^ of 
material within ~ 3 pc of the Galactic center. Since such a cluster would have experienced core collapse on 
about the same time scale, the most massive stars had already segregated to the cluster core. The deposited 
(core) material was therefore rich in massive stars. These findings are contrary to the results reported by 
Kim et al. (2002). 

More detailed studies are underway to qualify and quantify these statements (Portegies Zwart, 
McMillan & Gerhard 2003). Preliminary results indicate that the inspiral times derived here are in good 
agreement with A^-body calculations using the GRAPE-6 special-purpose computer (Makino et al 2002), 
with the same description of dynamical friction as presented here. More sophisticated calibration of the 
dynamical friction parameters themselves, obtained by modeling the Galactic background as individual 
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stars, will be the subject of a future paper. 
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Advanced Study and Tokyo University for their hospitality and the use of their GRAPE-4 and GRAPE-6 
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by the Space Telescope Science Institute and by NASA ATP grants NAG5-6964 and NAG5-9264, by the 
Royal Netherlands Academy of Sciences (KNAW) and the Netherlands Research School for Astronomy 
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Appendix A 

The argument X = Vc/V^cr in the dynamical friction relation (Eq. 7) of Sec. 2.1 may be evaluated as 
follows for inspiral through a sequence of nearly circular orbits. 

Following Binney & Tremaine (1987, Eq. 4-30), we write the equation of dynamical equilibrium (the 
radial Jeans equation) for stars near the Galactic center as 



where we assume an isotropic velocity distribution. In the power-law region, M oc R" (Eq. 1), we further 
assume that oc Vc"^. It then follows that a^p ~ so 



Substitution in Eq. (37) then yields X = v'2 — a. We note that, as a ^ 1, this reduces to the correct 
expression for an isothermal sphere (see Binney & Tremaine, p. 230). 



This preprint was prepared with the AAS lAT^jX macros v4.0. 
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